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ABSTRACT 



1 INTRODUCTION 



We assess the strengths and weaknesses of several likelihood formalisms, including the 
£h ■ XF a s t e r likelihood. We compare the performance of the XFaster likelihood to that of the 

ryj | Offset Lognormal Bandpower likelihood on simulated data for the Planck satellite. Parameters 

estimated with these two likelihoods are in good agreement. The advantages of the XFaster 
likelihood can therefore be realized without compromising performance. 
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. The temperature and polarization anisotropies of the Cosmic Microwave Background (CMB) contain a wealth of cosmological information. 
In extracting this information from measurements of the CMB, the likelihood function C(Ci) = P(data|C^ ), where Ci is the theoretical 
spectrum for some cosmological model, plays an important role. For Gaussian fluctuations, C(Ce) is given by a Multivariate Gaussian of 
the observed data (section |2j- For low resolution data (£ < 10 2 , defined more precisely later) in the usual spherical harmonic expansion, 
it is computationally feasible to evaluate this directly. For high resolution data, it is not, and computationally tractable approximations are 
' required. Fortunately, the large number of independent samples of the universe at high resolution, coupled with a much more nearly diagonal 
covariance matrix, mean that approximations exist that are accurate as well as fast. The subject of this paper is the performance of one such 
high I method, the "XFaster likelihood." 

A successful high t method must account correctly for correlations induced in the angular power spectra by partial sky coverage, non- 
uniform noise, and non-zero beamwidths, as well as the temperature and polarization cross power. A n umber of approaches for hi gh-l have 
been proposed. Fo r temperatu re alone, these include the Gauss ian, Offset Lognormal, Equal Variance l Bond. Jaffe. & Knoxll2000l) . WMAP 



( Verd e et al . 2003), and SCR ISmith. Challinor. & Rochj2006h ike lihood a pproxima tions. For temperature and polarizati on togethe r, these 



include the Offset Lognormal Bandpower, Hamimeche and Lewis l Hamimech e & Lewis! 20081) . and XFaster l Rocha et al. 20091) likeli- 
hoods. 

This paper is organized as follows. In Section [2] we describe the the Multivariate Gaussian likelihood both in pixel and harmonic 
space, give an overview of existing high £ approximations (Section [2. 2t and their limitations. An account of their performances applied to 
simulated Planck data has been given in jRocha et al . 2009). We furt her test the performance o f XFaster likelihood implemented in a new, 



modified version of the publicly available software CosmoMC coda 1 ! ( ILewis and Bridle 1 2002)) by applying it to estimation of cosmological 



parameters, and compare it to the Offset Lognormal Bandpower likelihood. We further test our approach to tackle the asymmetry of the beams 
by comparing parameters estimated for Planck simulated data convolved with a symmetric and a asymmetric beam using both likelihood 
approaches. 



1 http://cosmologist.info/cosmomc/ 
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2 LIKELIHOOD FOR A GAUSSIAN SKY 

Pixel temperature fluctuations T(n) (Stokes J) on the celestial sphere can be expanded in terms of spherical harmonics Yi m as 

T(n) = J2 a tmY em (n), (1) 

fm 

with coefficients a lm . Polarization fluctuations (Stokes Q and U) can be expanded in spin-2 spherical harmonics, 2 Ye m , with U and B (grad- 
or curl-type) polarization coefficients 

(Q ± ^)(n) = ^2(af m ± iaf m ) ±2 y fra (n). (2) 

em 

If the CMB is a Gaussian isotropic field, then the probability of a measurement of the sky given a model is described by a Multivariate 
Gaussian of the observed data: 

^lp)= 27rJV/2 | C|1/2 exp(-ldC-^) (3) 

where C is the covariance of the data d, and p is the set of model parameters. C(p) = S(p) + N, where S is the sky signal and N is the 
noise. Since measurements of the sky are pixelated, the above likelihood is often called the "pixel-based likelihood" when estimated in the 
pixel domain. It can be evaluated directly if the number of pixels is not too large, but for a full-sky experiment such as Planck with 5' pixels 
it is impossible with current computers. 



2.1 Exact likelihood in harmonic space 

If the CMB is Gaussian, its statistical properties are represented fully by the underlying power spectrum Ce- The multipole harmonic coeffi- 
cients, af m (where X is T, E, or B), on different scales are independent of one another, and we can write 

(af m ) afi m A = 5ai&mm'Cf x . (4) 



The af m are complex. Under the assumption of Gaussianity, their real and imaginary parts are independent and Gaussian distributed, hence 
their phases are random. Since the CMB is real, they must satisfy af^ — (— l) m af_ m . This means that af is real. 
For m = we have 

P(am\Ce)da m = ( 27r y V 2| C |i/ 2 ex P {-^oC^Vo} ■ (5) 
For m / 0, we have for the real part of the af m 

P(Vt{a em }\C e )dM{a em } = ^pps exp {-^a^jC^a^}} , (6) 
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and 

c t = i cr cr o i , (8) 

Similarly for the imaginary part. 

Combining together all the values of m, we find, for a particular £: 

P(C e \C e ) oc\C e \^\C e \-^ cxp(-^±iTr(cW) I. (9) 
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m=-e 

and the normalisation is independent of Ce and Ce- In data analysis, the measured power spectrum Ce is a fixed quantity, hence the 
dependence of the likelihood on this value is generally dropped. In this case, up to a constant, we can write the log-likelihood as 

- 2 In P {Ce | Ce) = (2£ + 1 ) (in | Ce. \ + Tr {Ce CJ 1 ) ) , (11) 

i.e., the inverse Wishart distribution. If we consider only one measurement, e.g., one T-mode or B-mode, we can write the likelihood function 
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( Bon d. Jaffe. & Knoxllioooh 



-21nP(C,|C,) = (2l + l) h n (g)+gj , (12) 

i.e., the inverse Gamma distribution, where Cj is the theoretical value of (or C e ) and CV is the measured value. 

We can write an exact expression for the likelihood function for our measured power spectrum Ce given the true underlying power 
spectrum Ce , which is a function of cosmological parameters. Since this likelihood is usually considered in the context of data analysis, it is 
common to regard the measured Ce as fixed quantities, and to write the likelihood as 

m P(C|C7) = ^-^±^ L (g) +gj + constant, (13) 

where the constant depends on Ce- For fixed Ce, this function peaks at Ce = Ce- However if we wish to consider the likelihood as a function 
of Ce then it is necessary to include the -dependence of the likelihood, in which case it should be written as 

In P(C\C) = Y (2£ ~ 1) In Ce - ^±3. ( m Cl + ^ j + constant. (14) 
- 2 2 y CeJ 

For a fixed underlying power spectrum Ce, this function peaks at Ce = ( 21+1 ) Ce- 

This is adequate for a full-sky experiment with an infinitely narrow beam and no instrumental noise. For a partial or "cut" sky it is 
necessary to account for the correlations between the Ce that are introduced. In addition, real experiments always have non-uniform noise, 
and must estimate bandpowers (Cs) rather than individual Ce- We need to find an appropriate likelihood function that includes the correct 
correlations and accounts properly for noise. 



2.2 Approximating the likelihood at high-^ 

Although the signal and the noise can be assumed Gaussian, the distribution of the band powers is non-Gaussian. This is the so called 
'cosmic bias' — the distribution is skewed towards higher Ci values. This effect is most noticeable at low i where cosmic variance of the 
signal dominates the error bars. As we will see, the full likelihood includes the cosmic bias whereas the Gaussian approximation of the Fisher 
matrix does not. The Joint likelihood for temperature and polarization carries an extra complication in that one has to find an approximation 
that properly accounts for the temperature and polarization cross power. We will start by considering the current approximations derived for 
temperature alone, followed by an account on existing attempts to extend it to a joint likelihood for temperature and polarization. 



2.2.1 Temperature only 

Gaussian Likelihood — The first level of approximation is to use a likelihood that is Gaussian in the Ce 1 Bond. Jaffe. & Knoxll2000l) . i.e., 

P(C\C) <x exp - CfS^iC - C) j , (15) 

where C is a vector of Ce val ues (and similarly C) and S -1 is the inverse signal covar i ance m atrix. However this likelihood function is well- 
known to be biased, (see e.g., Bond. Jaffe. & Knox 1 200C ); Smith. Challinor. & Rochal 1 20061) ). This Gaussian likelihood can be implemented 



in two ways. The version discussed in 



Bond, Jaffe, & Knox 



(2000) 



considers the case where the signal covariance matrix is derived from the 



measured power spectrum, rather than the theoretical power spectrum. This results in an overestimation of the errors if the measured power 
spectrum has fluctuated upwards, and hence upward fluctuations are given less wei ght, leading t o an o verall downward bias. The covariance 
matrix can instead be computed from the theoretical power spectrum, as used in 
amplitude. 



Verde et al 



(2003), leading to an overestimation of the 



Offset Lognormal Likelihood- 

2000), given by 



-A better approximation to the likelihood is the Offset Lognormal approximation ( IBond. Jaffe. & Knox 



i>L.xi.C\C) -v <-xp { ~(£ 



zY M(z 



where ze = ln(Ce + xe) and the matrix M is related to the inverse covariance matrix by 

Mur = {Ce + xe)See}{C'e + x'e) 



(16) 



(17) 



(The offset factors xe are simply a function of the noise and beam of the experiment.) This likelihood function is still slightly biased, but in 
the opposite direction to that of the Gaussian likelihood. 

To some extent the Offset Lognormal approximation addresses this issue, in that the transformation of variables from Ce to Ze = 
\n(C'e + xe) has a constant curvature matrix. Hence the uncertainties on the Ze do not depend on the ^.estimated, avoiding the cosmic bias. 
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To proceed, assume a normal distribution in this new variable Zt instead of Ce ■ The likelihood is now closer to the exact one. Nevertheless, 
there is still a small bias opposite to that introduced by the Gaussian approximation in Ce ■ 

WMAP Likelihood — Taking advantage of the fact that the bias for the Offset Lognormal likelihoo d is opposite to th at introduced by the 
Gaussian likelihood, the WMAP team defined a likelihood that is a weighted combination of the two ( Verde et al ]l2003h : 



lnP W MAp(C|C) = i In P G au SS (C I C) + |lnPL N (C|C) (18) 

This likelihood is a significantly better approximation for the case of a Gaussian CMB. However, the fact that a likelihood function is accurate 
in the absence of correlations (when the probability of a measured power spectrum is purely a function of the input power spectrum rather 
than having any additional dependence on the cosmology) is not a guarantee that it will perform well when applied to a non-Gaussian sky. 
For instance. lsmith. Challinor. & Rocha 1 20061) have shown that the WMAP likelihood gives significant biases in the dark energy parameter, 



w, when considering the lensed B-mode power spectrum. 

Equal Variance Likelihood — The equal variance likelihood proposed by Bond. Jaffe, & Knox 1 200d) is given by: 

lnL = -iG[e- ( -*> - (1 - (z - z))] (19) 



with 



and 



z = ]nlq b +q£ r ) (20) 



The noise offset is estimated using the equation of the maximum likelihood solution for the qb, replacing the observed map with the 
average of the noise Monte Carlo simulation power spectra < N, 



SCR Likelihood— ISmith. Challinor. & RochU (2006) developed a new likelihood to tackle the non-Gaussianity of the lensed sky for studies 



of the B-mode power spectrum; however, this likelihood has not been extended to the joint probability distribution for all modes. By 
considering the curvature (with respect to the measured Ce) of the exact log-likelihood expression for a Gaussian sky (equation 1 14b at its 
peak, and also the third derivative, a new likelihood can be derived which is Gaussian in xe = (Ce)^ , where both the second and third 
derivatives with respect to xe take the correct values at the peak of the likelihood. 

The SCR likelihood approximates the normalised distribution P(Ce\0) as Gaussian in some function of the Ce, and takes the form 



\nP{Ce\0) n]nA-\^M^(xe- Lie){x v - m), 



(22) 



where 

xe = C 1 / 3 (23) 
21 - 1 ^ 1/3 



and 

Here See' is the covariance matrix of the measured Ce at parameters 6. The normalisation is 

A' 1 oc (detM«0 1/2 II^' < 26 > 

e 

which can be approximated by A oc Yli 1/Ce. 

Applying this expression to a G aussian simulation gives results almost indistinguishable from the exact likelihood expression, as shown 
in Smith. Challinor. & Rochal 1 20061) . Ignoring the (21 — 1)/ (21 + 1) factors, however, is a bad approximation at low £, since it ignores the 



fact that, for a fixed Ce, the peak of the likelihood is slightly below Ce - This ends up translating to an underestimation of the theoretical power 
sp ectrum. The performance of the G aussian, WMAP, and two versions of the SCR likelihoods on full-sky lensed simulations was compared 
in lSmith. Challinor. & Rochal |2006). The WMAP likelihood function gives very different results to the new SCR likelihood function, and 



shows a significant bias in the dark energy parameter, w. 
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2.2.2 Temperature + polarization 



We cannot merely extend the above approximations to build a joint likelihood for temperature and polarization. For instance, assume that 
we approximate the likelihood for TE as a Gaussian. The mode and variance of the Gaussian distribution for TE depend on TT and EE. 
Hence it is not enough to consider the joint likelihood as a product of independent TT, EE, and TE likelihoods. Indeed the trick is to find 
a way of coupling these components reliably. Furthermore, given that the TE power spectrum is at times negative, we cannot build a Joint 
likelihood as a product of independent Offset Lognormal likelihoods. As a quick fix, one could try the following: 

Offset Lognormal Bandpower Likelihood — This likelihood is a joint likelihood for temperature and polarization built as a Gaussian for TE 
and Offset Lognormal for TT, EE, and BB. However, this approximation does not properly take into account temperature and polarization 
cross power. 

The following two likelihoods, the Hamimeche and Lewis and XFaster likelihoods, do attempt to take temperature and polarization 
cross power into account properly. 

Hamimeche and Lewis (HL) Likelihood — The Hamimeche & Lewis i 200 8h likelihood generalizes the full-sky exact likelihood given by 
Equation [TT] to the cut sky by considering a quadratic expression of the form 



\nL(Ct\C e ) = [ 9 {D tiU f = ^±l T r [g(D« 



(27) 



where g(x) — sgn(a; — l)y / 2(x — In x — 1) and [g(Df)] i;; = g(De.u)8ij . This approximation involves a fiducial model so that the 
covariance can be pre-computed. It is assumed that the matrix of estimators Ce is positive/definite, although this assumption may break 
down for some estimators at low-£. 

XFaster Likelihood — The XFaster likelihood, introduced in lContaldi et al. | J2009I) and lRochaetal] J2OO9I) . takes the following form for 
temperature alone: 



r^obs 



(Ce + (tie)) 



+ ln 



(pi + (fit)) 



(28) 



where a tilde connotes a quantity estimated on the cut-sky. The power spectrum is parameterized through a set of deviations qe from a 
template full-sky spectrum C x 



(S) 



(29) 



where Kee 1 is the coupling matrix due to the cut sky observations, Be expresses the effect of a finite beam, and Ft is a transfer or filter 
function accounting for the effect of pre-filtering the data in both time and spatial domains. 
Extending to polarization, we have 



where the matrix C is block diagonal: C 
matrix 



lni = -2^3(2^ + 1) (Tr(Df 
diag(Y)e n 



D, 



He + CN t 
ii • • • , D«„ 



^e 



cf 



ci 

cf 



+ In Df + (N 



(30) 



with each multipole's covariance given by the 3x3 



(31) 



This likelihood follows intuitively from the full-sky exact likelihood, the Inverse Wishart distribution, as given by Equation [TT1 



2.3 The likelihood at low- 



The pixel-based Multivariate Gaussian likeliho od given by Equation [3] can be co mputed up to £ ~ 30, 4 0, and is adequate for comparison 



with Xf aster as shown in jRocha et alj2009h . (Faster methods (see summary by 



Ashdown et al 



1 20 id )) have been developed, but 



are not 



necessary here.) We use an implementation known as Bf like, described as follows. 

A CMB map can be written as an ordered vector d = (T 1 , T 2 , T, lT , Qj 1 , Qj 2 , ■■•Qn P , Uh , Uj 2 1 —U np ), comprising all pixels 
with valid observations. In general tit 7^ np and the sets of indexes of temperature and polarization measurements will be different. 
Assuming that both CMB and noise fluctuations in each pixel are Gaussian-distributed with zero mean, the likelihood for d has the functional 
form given in equation (1), where the covariance matrix has a block structure: 

( < TT >(n T xn T ) 



C = 



< TQ > (n 



(np X tit) 



(tit Xnp) 

< QQ >( npXnp ) 



<TU>, n 



(np Xnp) 

<QU> 

(n p Xnp) 

<UU > 

(np Xnp) , 



(32) 



6 G. Rocha et al. 



Correlations between, e.g., temperature measurements in pixels i\ and 12 can be written as: 

(T n T t2 ) = £ ^±lc e P e (9 ni2 ) +N, lla , (33) 

Pi(x) are the ordinary Legendre polynomials, and ft;^ is the angle between the centers of pixels ii and 12. Notice that the {Ci} include 
the contribution of the beam and pixel window, i.e., {Ci} = {C| h }fc|iu|, where {C^ h } is the theory power spectrum and bi and wt are 
the harmonic transform of the beam and window functions respectively. For uncorrelated noise, Nij = nfSij , where rij is the rms noise in 
pixel i . In general, N is a dense matrix. Similar expressions hold for correlations involving Q and U (see e.g. jTegmark & de Oliveira-Costal 
(2001)). The choice of ^ max in Equation [33] depends on several factors, including the smoothing beam, the signal-to-noise ratio, and the 
pixelization scheme. 



3 COMPARING LIKELIHOODS 
3.1 Simulations 



To compare th e performance of XFaster with other likel ihood functions, we use simulations developed within the Planck CTP working 
group. Planck 1 Planck Blue Book 2005 : Tauber et al. 2009) is a full-sky expe riment covering frequencies from 30 to 85 7 GHz with beams 



ranging in size from 33' to 5'. A full description of the simulations is given in lAshdown et al. 1 20ld) : 



Rocha et al 



(2009). 



Practical consid- 
erations of computational resources having to do with the size of the time-ordered data (TOD), the number of pixels in the maps, and the 
number of multipoles that had to be calculated, dictated the choice of the 70 GHz channel for the simulations. Higher frequency channels have 
higher angular resolution and sensitivity, and will extend to smaller angular scales with reduced error bars. Recent increases in computational 
capability make it pos sible now to generat e thousands of Monte C arlo simulations at the higher frequencies as well. Results will be presented 
in a future publication IRocha et al. (2010): Ashdo wn et al. 

The 70 GHz simulations used here include the CMB, realistic detector noise, and noise induced by temperature fluctuations of the 20-K 
hydrogen sorption cooler. To mimic the sensitivity of combination of channels, as would be used for separation of foregrounds with real data, 
the white noise level was taken to be lower than that expected for the 70 GHz channel alone. The white noise per sample was 2025.8 fiK and 
the 1// noise power spectrum had knee frequency 0.05 Hz and slope —1.7. The input sky signal used to generate the "observed map" was 
the CMB map derived from the Planck CMB reference sky available irjf] which uses the WMAP 1-year aim up to i = 70 to generate the 

CMB. In other words, the large scale structure of the observed map is a WMAP-constrain ed realization. 

The TOD were generated using modules of the Planck simulator pip eline, Levels feeinecke et al.||2005l). Maps were made from th e 
simulated time-ordered data (TOD) using the destriping code Springtide dPoutanenll2005l : lAshdown et al.ll2007a lbi l2009l : lAshdownll20 09b). 
Where a sky cut was applied, it was made at the boundary where the total intensity of the diffuse foregrounds and point sources exceeded 
twice the CMB sigma. Pixels missing due to the scanning strategy were masked. The beams of the detectors have FWHMs of 13'— 14', so the 
maps were made with N s - 1( ic ~ 1024, corresponding to a pixel size of 3.'4. Two sets of maps were provided, one 12-detector map to be used 
in the auto-spectrum mode, and three 4-detector maps to be used in the cross-spectrum mode. 

Two cases were considered. The first, called Phase 2a for historical reasons, assumed symmetric Gaussian beams with FWHM of 14'. 
The second, Phase 2b, assumed elliptical Gaussians fit to the central parts of realistic beams calculated by a full diffraction code for the 
Planck optical system. 

One hundred Monte Carlo signal simulations were generated from the best fit WMAP + CBI + ACBAR ACDM power spectrunjfl with 
BB mode power set to zero. For the symmetric beam case the signal could be simulated in the map domain, while noise was generated in the 
time domain. For the asymmetric beam case, both signal and noise simulations were generated in the time domain. 

The "low-^ dataset" of the Phase2 simulations was generated directly at A^dc = 16. The procedure adopted ensured consistency 
between the low- and high-f datasets used to test the XFaster power spectrum and likelihood estimator; however, the low-£ dataset lacks 
the artifacts connected t o smoothing and de gradation of higher resolution maps that will be present in the final Planck maps. 

As pointed out in IRocha et alj ( l2009h . since the large scale structure of the observed map is derived from real observations, i.e., a 
WMAP constrained realization, it is not necessarily consistent with the best-fit spectrum at low multipoles. This discrepancy is evident 
when comparing the cosmological parameters estimated with XFaster power spectrum and likelihood and the Offset Lognormal likelihood 
with the theoretical best fit parameters. The Monte Carlo simulations, however, are realizations of the first year WMAP+CBI+ACBAR best 
fit ACDM power spectrum for Phases 2a and 2b, so such discrepancy is no longer present. Parameters estimated from these Monte Carlo 
simulations maps are shown to be close to the WMAP best fit parameters. 



3.2 Comparisons 

Rocha etal.|j2009h (Figs. 14 and 15) showed that for TT and TE, all of the likelihood approximations described above except the Gaussian 



2 http://www.sissa.it/ plancWreference-sky/CMB/alms/alm-cmb-reference-template-microKthermodynamic-nside2048.fits 

3 available in http://lambda.gsfc.nasa.gov/product/map/drl/lcdm.cfm 
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likelihood converge to the same form at multipoles above 10 or 20. For EE, however, the Gaussian and the Lognormal likelihood differs 
from the rest up to high I ~ lO. lRochaetai] 12009) (Fig. 17) also showed by compar ng XFaster with the pixel-based likelihood code 
(BFlike) that a transition from low-f to high-£ codes was appropriate for Planck in the range ftrans = 30-40 

Here we assess the performance of the XFaster likelihood by comparing cosmological parameters obtained with the XFaster 
and the Offset Lognormal Bandpower (i.e., Offset Lognormal likelihood for TT, EE, BB, and Gaussian for TE) likelihoods. XFaste r 
computes the likelihood of a model passed to it by a modified version of the publicly available CosmoMC cod^f] ( TLewis and Bridle i 2002h )) 
for cosmological parameter Markov Chain Monte Carlo estimations. There is no need for window functions or the band power spectrum 
itself. The inputs are the raw pseudo-CV of the observations plus the kernel and transfer function required by XFaster to relate the cut-sky 
pseudo-CV to the full-sky CV. 

For the Offset Lognormal Bandpower likelihood, window functions are required that properly accoun t for band power spectrum cor- 



relations. We used two different window functions, top hat (box) and Fisher-w eighted jRocha et al. 200S), writt en Tbb- For the sak e of 



completeness we describe here how we derive Tbb- To construct Tbb ( following iBond. Jaffe. & Knox (2000) and lRochaetai] l2009h ) 
define a logarithmic integral, 

£■ 



t{t + 1) 



(34) 



which is used to calculate the expectation values for the deviations qf,(when a shape model, Cf is considered), or bandpowers Cb (when Cf 
is assumed to be flat). 



(Q b ) = 



T\W e b C el 



where W\ is the band power window function, with C^ s ' = £(£ + 1)C I 



(C b ) = 



T [WjCt] 

i \wf\ 



/2tt. We define normalized window functions 

= 1, 



using the fact that 



to obtain 



Extending to polarization: 



((Cf s - N e )) Ct 



4tt 



(2£+l 



■E^' 1 E^' + 1 ) 



^b'i 



(S) 



(cf, + (N e ,)y 



-Ktt'FiBi . 



(21 + 1) 



E E + ^ Tr l w »'t> K * 



where Wi 



D 



-l as -ft-i 



and K.e gives the cut-sky response to the individual full-sky multipoles, 



/ K ei Ff T B 



,K„ t FT E B, 



\ 



c Ki> iFj B ' B\ 



-K t 
(- 



x KinF e 

pEE R 2 , 

-Kiii —- Kff 



Ki>iFf B Bf 
Ff B Bp 



Ki'iFf 3 Bj 



? EB 



-K e , e F f BB B] + 



K tu ) F L 



(35) 



(36) 



(37) 



(38) 



(39) 



(40) 



B 



The parameters calculated are: baryonic, cold dark matter, and cosmological constant densities, ojb = Qbh 2 



— Q c h 2 , and uja 



ft Ah 2 , respectively; the ratio of the sound horizon to the angular diameter distance at decoupling, 8 a ; the spectral index of the initial 
fluctuation spectrum, n B ; the overall normalization of the spe ctrum 1ok[10 10 A1 a t k — 0.05 Mpc -1 (A s ); the Hubble constant Ho; and 
the reionization redshift z re . We treat r in two different ways. iRocha et al.U2009l) showed that "high €' codes could be used to determine 
parameters from the "observed map" quite well if r was fixed in the fit to the value of the input model. This works because r is constrained 
primarily by data at £ < 30. We indica te when r is held fixed. 

FigureU from Roc ha et al. i 2009h . shows one-dimensional marginalised parameter distributions from Xf aster, for three cases: 1) the 
observed map; the ensemble average of Monte Carlo simulations; and the observed map, but holding r fixed at its input value and using the 
Xf aster likelihood only for £ > 30. The input parameters are recovered quite well from the ensemble average. The red lines show that if 
t is fixed, high-^ codes such as XFaster can ignore the low multipoles that they are not designed to calculate, and still recover the other 
input parameters quite well. 

Figure|2]shows parameter distributions from the Offset Lognormal Bandpower likelihood for a Fisher-weighted (Tbb) window function, 
for both the observed map and for the ensemble average of Monte Carlo simulations. The input parameters are recovered from the ensemble 
average simulated data, but not from the observed map, particularly A s . This is a not a surprise. As described in § 3.1, the observed map 
is fixed by WMAP for 2 ^ I ^ 70. It is therefore affected by low-^ anomalies arising from any cause, including the details of WMAP 



4 http://cosmologist.info/cosmomc/ 
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Table 1. Par ameter estimates an d uncertainties from the Offset Lo gnormal Ba ndpower and XFaster likelihoods for Planck simulations at 70 GHz, compared 
to WMAP (Dunklev et al. 2009) and Fisher uncertainties (Rocha et al. (2004)). Estimates are for the ensemble average of 100 Monte Carlo simulations. Last 
column displays input parameter values ± Fisher uncertainties for reference. 



Parameter WMAP Offset Lognormal Bandpower XFaster Fisher 



T 


0.087 ±0.017 


0.099 ±0.015 


iinti+ 00643 
u.liuo_ 00771 


0.1103 ±0.004 (4%) 


rt s 


0.963 ±0.015 


0.965 ± 0.008 


o 6 2l+ - 01130 
u - auz -0.01170 


0.9582 ± 0.004 (0.4%) 


Ub 


0.02273 ±0.015 


0.0229 ± 0.0003 


0.0225 ± 0.00042 


0.02238 ± 0.00018 (0.8%) 



processing. The WMAP best-fit parameters, on the other hand, are obtained with heavy marginalization of the low-f points by foregrounds, 
and are therefore little affected by the low-<? anomalies. Since the Monte Carlo simulations are realizations of the WMAP best-fit model 
parameters, we expect no systematic bias from the ensemble of simulations, but significant offsets from the parameters derived from the 
observed map, as confirmed. 

Figure [3] shows the effect of fixing r in the case of the Offset Lognormal Bandpower likelihood for a top-hat window function for the 
observed map. Green dashed lines are computed from all multipoles, with r free to vary, while for the blue solid lines use only I > 30, with 
r fixed at the input value. As with XFaster in FigureQ] fixing r and ignoring low multipoles gives good results with the other parameters. 

Figures[4]and[5]show parameter distributions from the Offset Lognormal Bandpower likelihood in the symmetric and asymmetric beam 
cases, for the observed power spectrum and the ensemble average of the Monte Carlo simulations, respectively. Parameters for the two cases 
are roughly cons istent with each other. They can be compared to the equivalent distributions from XFaster likelihood in Figure 21 of 
Roch a et al. I J2009I) . 

Investigating the plot for the average mode, for the case XFaster likelihood we see deviations of the order of cr/2 for Q c h 2 , as, n s 
and Ho. There is an obvious degeneracy between a$ and n s . For the observed case these deviations are noticeable mostly in A s and a$. 
Once again these parameters are degenerate. 

The overall agreement in the parameter constraints from both symmetric and asymmetric beams is quite impressive. This reflects the 
adequacy of our procedure when dealing with beam asymmetries, although there is still a slight bias for the asymmetric beam case. This bias 
is consistent with the small differences in the estimated parameters obtained with the Offset Lognormal Bandpower likelihood, in particular 
for parameters such as n B , as and log[10 10 A s ]. Although the power spectra look consistent, the parameter estimation shows departures of 
the order of a/ 2 for some of the parameters. 

XFaster assumes that the noise is white (uncorrelated), i.e., that the noise covariance matrix is diagonal. Also, the XFaster likelihood 
is estimated multipole by multipole, unlike the Offset Lognormal Bandpower likelihood. Hence a proper estimation of the transfer function 
(filter function) requires a large number of Monte Carlo simulations to reduce the correlations between multipoles introduced by, e.g., sky 
cuts. These simulations include both correlated noise and a sky cut. Hence for the small set of 100 Monte Carlos used here we should expect 
a larger deviation when XFaster likelihood is employed. 

Figures[6]|7] and[8]compare the effect of top-hat and Fisher (J-'bb) window functions on parameters estimated with the Offset Lognormal 
Bandpower likelihood. Fisher-weighted window functions account for the band power spectrum correlations. It is clear that most parameters 
imp rove with Tbb wind ows, while uncertainties for all but r are unaffected. The 95% upper limit on r is higher by 15% for J-'bb windows. 

J2009I) ameter constraints obtained with XFaster likelihood for the symmetric beam case. TableQ]is the 
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same for the Offset Lognormal Bandpower likelihood and the XFaster likelihood for the ensemble average power spectrum of the Monte 
Carlo simulations. O ur aim is solely to check how the parameter uncertainties for Planck from both likelihoods compare to those of WMAP 
1 Dunklev et"al]|2 009) and to our Fisher predictions. The uncertainties are 2-3 times better than those for WMAP except for r estimated with 
the Offset Lognormal Bandpower likelihoo d, for which the uncertainty on r is 0.015 compared to ~0.007 for XFaster, 0.017 for WMAP, 
and 0.004 for our Fisher predictions dRocha et alJ II2004I) . Table[T|l. 



4 CONCLUSIONS 

Parameters estimated with the XFaster and Offset Lognormal Bandpower likelihoods agree well. As the XFaster likelihood is estimated 
for individual multipoles, a large number of Monte Carlo simulations is required for accurate estimates of low-^ correlations. If only a small 
number of Monte Carlo simulations (such as the 100 used in this study), binning of the band power spectrum estimated with XFaster 
used along with the Offset Lognormal Bandpower likelihood helps to partially correct these correlations. For a large number of Monte Carlo 
simulations this is unnecessary. There XFaster likelihood, however, has at least three advantages. First, the Offset Lognormal likelihood 
does not properly take into account the temperature-polarization cross power. This is likely to become evident with a larger number of 
simulations, or at lower nois e level s, such as anticipated with the Planck HFI 143 GHz channel. We are investigating this further; results will 
be presented in lRocha et al. (2010). Second, the Offset Lognormal Bandpower likelihood requires calculation of a window function. Third, 
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Figure 1. Marginalised parameter distributions from Xfaster, for the observed map (blue dashed lines), the ensemble average of Monte Carlo simulations 
(black solid lines), and for the observed map holding r fixed at its input value and using the Xfaster likelihood only for i > 30 (solid red lines). Input 
parameter values are marked by vertical lines. The input parameters are recovered quite well from the ensemble av erage (see text), and also for the observed 
map when r is fixed and low multipoles are not included in the parameter fits. Figure reproduced from Rocha et al. (2009) 



the XFaster likelihood can go straight from maps to parameters (via its raw pseudo-Cf), bypassing the band power spectrum estimation 
step. These advantages make XFaster a adequate procedure to estimate cosmological parameters from Planck data in the high multipole 
regime. As a bonus, XFaster performs reasonably well for moderately low multipoles as well. Although hybridization with a likelihood 
code able to handle fully the challenges of multipoles less than, say, 40, will be necessary for the best estimates parameters, XFaster could 
be used alone where accuracy can be traded for speed. 
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Figure 2. Marginalised parameter distributions from the Offset Lognormal Bandpower likelihood (Offset Lognormal for TT, EE, BB, Gaussian for TE), 
for a Fisher- weighted (F^b) window function. Blue solid lines are from the observed map. Red dashed lines are from the ensemble average of Monte Carlo 
simulations. Input parameter values are marked by vertical lines. 
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Figure 6. Comparison of different window functions with the Offset Lognormal Bandpower likelihood. Parameter distributions are shown for the observed 
map for top-hat (green dashed lines) and a Fisher (Tbb) (dark blue solid lines) window functions, for the asymmetric beam case. Input model parameter values 
are marked by vertical black lines. Most parameters improve with J 7 ;,;, -windows, while uncertainties are unaffected except in the case of r. The 95% upper 
limit on r is higher by 15% when using J 7 ^ windows. 
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Figure 7. Parameter constraints from the Offset Lognormal Bandpower likelihood with a top-hat window function, computed for the observed power spectrum 
(black dashed lines) and for the ensemble average of 100 Monte Carlo simulations (blue solid lines) 
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Figure 8. Same as Figure|7]for a Fisher (Fbb) window function, computed for the observed power spectrum (blue solid lines) and for the ensemble average of 
100 Monte Carlo simulations (red dashed lines). 



